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Abstract 

We prove that the ground state energy of an atom confined to two dimensions with 
' an infinitely heavy nucleus of charge Z > and iV quantum electrons of charge —1 

is E(N,Z) = -\Z 2 InZ + (E TF (X) + ±c H )Z 2 + o{Z 2 ) when Z ->■ oo and N/Z -> A, 
where E TF (X) is given by a Thomas- Fermi type variational problem and c H ~ —2.2339 
is an explicit constant. We also show that the radius of a two-dimensional neutral 
atom is unbounded when Z — > oo, which is contrary to the expected behavior of three- 
dimensional atoms. 
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1 Introduction 



We consider an atom confined to two dimensions. It has a fixed nucleus of charge Z > 
and N non-relativistic quantum electrons of charge —1. For simplicity we shall assume that 
electrons are spinless because the spin only complicates the notation and our coefficients in 
an obvious way. The system is described by the Hamiltonian 

= E (-^ ~ ifr) + E t^-\ 

i=l v 1 11 7 l<i<j<N 1 J 1 

acting on the antisymmetric space /\t =1 L 2 (R 2 ). Note that we are using the three-dimensional 
Coulomb potential to describe the confined atom. The ground state energy of the system is 
the bottom of the spectrum of H N Zl denoted by 

E(N, Z) = inf spec H N<Z = inf (ip, H N ^). 

W 1,2=1 

One possible approach to obtain the above Hamiltonian is to consider a three-dimensional 
atom confined to a thin layer IR 2 x (—a, a) in the limit a — > + (see [3], Section 3, for a detailed 
discussion on the hydrogen case). 

To the best of our knowledge, there is no existing result on the ground state energy and 
the ground states of the system, except for the case of hydrogen [221 [EH]. The purpose of 
this article is to give a rigorous analysis for large Z-atom asymptotics and our main results 
are the following theorems. 

Theorem 1 (Ground state energy). Fix A > 0. When Z — > oo and N/Z — » A one has 
E(N, Z) = - l -Z 2 lnZ+ (V F (A) + ^ Z 2 + o(Z 2 ) 

where E TF (X) is the Thomas-Fermi energy (defined in Section^) and c H = — 31n(2) —27^ + 
1 rs —2.2339 with 7^ f» 0.5772 being Euler's constant 0]/. In particular, A 1— > £7 TF (A) is 
strictly convex and decreasing on (0, 1] and E TF (X) = E TF (1) if A > 1. 

Remark. By considering the hydrogen semiclassics we conjecture that the next term of 
E(XZ,Z) is of order Z 3 ^ 2 . In contrast, the ground state energy in three dimensions be- 
haves as 

E(Z, Z) = c TF Z 7 ^ + c s Z 2 + c DS Z^ 3 + o(Z 5 / 3 ), 

where the leading (Thomas- Fermi [281 [6]) term was established in [15], the second (Scott 
|21j ) term was proved in [23], and the third (Dirac-Schwinger [2j [20]) term was shown in 

For three-dimensional atoms the leading term in the energy asymptotics of order Z 7 ^ 3 
may be understood entirely from semiclassics. The contribution to this term comes from 
the bulk of the electrons located mainly at a distance of order Z~ 1 / 3 from the nucleus. The 
term of order Z 2 , the Scott term, is a pure quantum correction coming from the essentially 
finitely many inner most electrons at a distance of order Z~ l from the nucleus. 

In the two-dimensional case the situation is more complicated. The leading term of 
order Z 2 ln(Z) is semiclassical and comes from the fact that the semiclassical integral is 
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logarithmically divergent, but has a natural cut-off at a distance of order Z~ x from the 
nucleus. The term of order Z 2 has two contributions. One part is semiclassical and comes 
essentially from electrons at distances of order 1 from the nucleus and another part, which 
corresponds to the three-dimensional Scott correction, coming from the essentially finitely 
many inner most electrons at a distance of order Z~ x from the nucleus. 

Thus the two-dimensional atom has two regions. The innermost region of size Z~ x 
contains a finite number of electrons and contributes with Z 2 to the total energy. The 
outer region from Z~ x to order 1 has a high density of electrons and can be understood 
semiclassically. It contributes to the energy with both Z 2 ln(Z) from the short distance 
divergence and with Z 2 from the bulk at distance 1. 

Theorem 2 (Extensivity of neutral atoms). Assume that N/Z — > 1 and ^n,z is a ground 
state of Hn z- Then, for any R > there exists Cr > such that 



/ p^ NZ {x)dx > C R Z + o(Z). 

J\x\>R 



Remark. If we define the radius Rz of a neutral atom (N = Z) by 

p^ zz {x)dx = 1 

\x\>R z 

then Theorem [2] implies that lim^oo R z = oo. In three dimensions, however, the radius is 
expected to be bounded independently of Z (see (221 [21]) • 

Our main tool to understand the ground state energy and the ground states is the 
Thomas-Fermi (TF) theory introduced in Section El In this theory, the Z-ground state 
scales as ZpJK z (x) and the absolute ground state pJ F (when N = Z) has unbounded sup- 
port. Roughly speaking, the extensivity of the TF ground state implies the extensivity of 
neutral atoms (in contrast, the three-dimensional TF Z-ground state scales as Z 2 p TF (Z 1 ' 3 x), 
i.e. its core shrinks as Z^ 1 ^ 3 ). 

The challenging point of the two-dimensional TF theory is that the TF potential V TF (x) 
is not in L 2 oc (IR 2 ) (it behaves like near the origin). Consequently, one cannot write the 
semiclassics of Tr [— h 2 A — V TF ] _ in the usual way because 

w// lh y - = J [v-wr* = -oo. 

This property complicates matters in the semiclassical approximation. In contrast, the three- 
dimensional semiclassical approximation leads to the behavior — (157r 2 /i 3 ) -1 J r3 \V TF (x)]+ 2 dx 

which is finite for the Coulomb singularity V TF (x) ~ G (R 3 ). 

We shall follow the strategy of proving the Scott's correction given by Solovej and Spitzer 
[2S] (see also [2"5]). that is to compare the semiclassics of TF-type potentials with hydrogen. 
More precisely, in the region close to the origin we shall compare directly with hydrogen, 
whereas in the exterior region we can employ the coherent state approach. We do not use 
the new coherent state approach introduced in [26], since the usual one [91 [27] is sufficient 
for our calculations. In fact, we prove the following semiclassical estimate for potentials with 
Coulomb singularities. 



3 



Theorem 3 (Semi classics for Coulomb singular potentials). Let V G L 2 OC (IR 2 \{0}) be a 
real-valued potential such that ln x \>nV + G L 2 (R 2 ) and 

\V(x) - nix]- 1 ] < C\x\- 9 for all \x\ < 5, 

where k>0,5>0,1>6>0 and C > are universal constants. Then, as h — > + , 

Tr[-h 2 A-V]_ = -(87T/1 2 )- 1 / {[V(x)} 2 + - k 2 ^ 1 - lf + ) dx 

Jr 2 

+K 2 (Ah 2 )~ l [hi{2K- l h 2 ) + c H ] + o(h- 2 ), 
where c H = -31n(2) - 2^ E + 1 w -2.2339 with ^ E ~ 0.5772 being Euler's constant f$. 

The article is organized as follows. In Section [2] we give a brief summary of the existing 
results concerning atoms confined to two dimensions. Section [3] contains basic information 
on the TF theory. The most technical part of the article is in Section HI where we show 
the semiclassics for the TF potential. The main theorems are proved in Section [51 Some 
technical proofs are deferred to the Appendix. 

2 Preliminaries 

2.1 Spectral Properties 

For completeness, we start by collecting some basic properties of the spectrum of H^,z, whose 
proofs can essentially be adapted from the usual three-dimensional case (see the Appendix). 

Theorem 4 (Spectrum of Hn,z)- Let H^,z be the operator defined above. 

(i) (HVZ Theorem) The essential spectrum of H^,z is 

ess spec Hn z = [E(N — 1,Z), oo). 

Consequently, for non-vanishing binding energy E(N—1)—E(N) =: e > 0, the operator 
Hn,z has (at least) one ground state. Moreover, in this case any ground state ^n,z of 
Hn,z has exponential decay as 

P^ N .z( x ) < C\x\^' 2 e' 2 ^ 2 ^ for \x\ large 

where the density p^ N z is defined as in Section \2.2\ 

(ii) (Zhislin's Theorem) If N < Z + 1 then the binding condition E(N) < E(N — 1) is 
satisfied, and hence H^,z has a ground state. 

(Hi) (Asymptotic neutrality) The largest number N = N C (Z) of electrons such that H^^z 
has a ground state is finite and satisfies lirnz-^oo N C (Z)/Z = 1. 

In particular, the spectrum of hydrogen (N = 1) is explicitly known [2H] (see also [TH] 
for a review). 
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Theorem 5 (Hydrogen spectrum). All negative eigenvalues of the operator — ~A — \x\ 1 in 
L 2 (R 2 ) are 



E n - 



Tr 



2(n + 1/2) 2 ' 
with multiplicity 2n + 1, where n — 0, 1, 2, ... 

The following consequence will be useful in our estimates. The proof can be found in the 
Appendix. 

Lemma 6 (Hydrogen semiclassics) . When fi — > + we have 

±A - {xl- 1 +fM = l - [ln<» - 3 ln(2) - 2 lE + 1] + o(l). (1) 

By scaling, for \i > fixed and h — > + , 

Tr [-/i 2 A - + /x] _ = (4/i 2 )^ 1 [ln(2/i 2 ) + ln(/i) + c H ] + o(/T 2 ), (2) 
where c H = -31n(2) - 2-f E + 1 w -2.2339 with ^ E ~ 0.5772 beinp £«/er's constant Jffl. 

2.2 Useful Inequalities 

For the readers' convenience, we recall some usual notations. We shall denote by L 2 (IR 2 ) 
the Hilbert space with the inner product (/,<?) = J R2 f(x)g(x)dx. An operator 7 on L 2 (IR 2 ) 
is called a (one-body) density matrix if < 7 < 1 and Tr(7) < 00. Its density is Pj(x) := 
j(x,x), where j(x,y) is the kernel of 7. More precisely, if 7 is written in the spectral 
decomposition 7 = J^i^i \ u i) tnen 7(^)3/) := Si Uui{x)Ui(y) and p 7 (x) := ^ ^^(a;)! 2 . 
For example, the density matrix 7^ of a (normalized) wave function ^ G Aili -^ 2 (^ 2 ) is 



7*(x, y) := N ^(x, x 2 , x N )^(y, x 2 , x N )dx 2 ...dx N , 

JR2(JV-1) 

which satisfies < 7$ < 1 and Tr(7^) = iV. Moreover, its density is 

Pm,(x) := p 7<p (a;) = N / |^(a;, x 2 , x N )\ 2 dx 2 ...dx N . 

JR2(JV-1) 

The following theorem regarding the spectrum of Schrodinger operators is important for 
our analysis (see e.g. [12] for a proof). The analogue in three dimensions was first proved 
by Lieb and Thirring [T7] . 

Theorem 7 (Lieb-Thirring inequalities). There exists a finite constant L 12 > such that 
for any real-valued potential V with V+ G L 2 (M. 2 ) one has 

Tr[— A - V\- > -Lx. 2 / V 2 (x)dx, (3) 

where a + := max{a, 0} and a_ := min{a, 0}. Hence Tr[— A — V]- is the sum of all negative 
eigenvalues of —A — V in L 2 (IR 2 ). 

Equivalently, there exists a finite constant K 2 > such that for any density matrix 7 one 

has 

Tr[-A 7 ] > K 2 ! p 2 Jx)dx. (4) 
Jr 2 
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Note that in general there is no upper bound on Tr l(_ OO)0 i(— A — V), the number of 
negative eigenvalues of —A — V, in term of J V? for any a > 0. However, we shall only need 
some localized versions of this bound. The proof of the following lemma can be found in the 
Appendix. The estimate in (ii) is useful to treat the Coulomb singularity in the region close 
to the origin (recall that ^ L 2 oc (IR 2 )). 

Lemma 8. Let V : M 2 — > R and let < <j){x) < 1 supported in a subset Q C M 2 with finite 
measure \Q\. Let < 7 < 1 be an operator on L 2 (M. 2 ) such that 

Tr[(-/i 2 A - V)(f)^(j)} < for some 1/2 > h > 0. 

(i) If V+ G L 2 oc (M. 2 ) then (jyycj) is trace class and there exists a universal constant C > 
(independent ofV, 7 and h) such that for any a G [0, 1], 

/ ^(^^CA-^iiy+ii^njini 1 -*. 

(ii) If V(x) < Co(|x| _1 + 1) then 070 is irace c/ass and i/iere exisis a constant C > 
dependent only on Cq (but independent ofV, (ft, Q, 7 and stzc/i that for any a G [0, 1] ; 

/ P ^(x)dx < ch- 4a (| in fc] + \n\) a in] 1 -*. 



We shall approximate the ground state energy E(N, Z) by one-body densities. For the 
lower bound, we need the following inequality to control the electron-electron repulsion 
energy. The three-dimensional analogue of this bound was first proved by Lieb [TU] and was 
then improved by Lieb and Oxford [13]. The two-dimensional version below was taken from 

Theorem 9 (Lieb-Oxford inequality). For any (normalized) wave function \I> G Ai=i -^ 2 (^ 2 ) 
it holds that 

[% E r4 — r* ) ^ D M~ c ^o [ pT , 

V l<i<j<N X 3\ ) J 

with Clo — 192(27r) 1//2 ; where the direct term D(p^) is defined as in Section HOI 

For the upper bound, we shall need the next result [TT] . 

Theorem 10 (Lieb's variational principle). For Z > 0, N G N and any density matrix 7 
with Tr(7) < iV ; one has 



E(N, Z) < Tr 



-A - Z\x\ 



1 



+ D{p~) - 



\l(x,y)\ : 
\x - y\ 



dxdy, 



where the direct term D{p 1 ) is defined as in Section HOI 
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2.3 Coulomb Potential 



Here we study the Coulomb potential / * | . | 1 of some function /. Associated to this 
potential is the Coulomb energy of two functions, 



DU,9)-All } -P^^V. 



2 J J \x — y\ 

R 2 xR 2 

That D(f,g) is well-defined at least in L 4 / 3 (IR 2 ) is due to the Hardy-Littlewood-Sobolev 
inequality (see [12], Theorem 4.3) 

D(\f\, \g\) < C HLS \\f\\ Li/3 \\g\\ Li/3 for all f,g G L 4 / 3 (R 2 ). 

Moreover, \x — y\~ l is a strictly positive kernel since the 2D Fourier transform of | . | 1 
is itself up to a constant (see [12] Theorem 5.9). Therefore, D(f) := D(f,f) is always 
nonnegative and (/, g) i— > D(f, g) is a positive inner product in L 4 / 3 (R 2 ). These observations 
allow us to formulate the following theorem. 

Theorem 11 (Coulomb norm). There exists Chls such that 

< £>(/) < Chls ||/||£v3 for all f G L 4/3 (M 2 )\{0}. 
Consequently, f i-> ^jD(f) is a norm in L 4 / 3 (M 2 ). 

In three dimensions, the Coulomb potential p * \ . | _1 of a radially symmetric function p 
is represented beautifully by Newton's Theorem (see [12J, Theorem 9.7). In two dimensions, 
however, we do not have such a representation since | . | _1 is not the fundamental solution 
to the two-dimensional Laplace operator. Therefore, the following bounds will be useful in 
our context and their proofs can be found in the Appendix. The lower bound is similar to 
Newton's Theorem in three dimensions, but the upper bounds are more involved. We do not 
claim that they are optimal but they are sufficient for our purposes. 

Lemma 12 (Coulomb potential bound). Assume that p is radially symmetric, < p{x) < 
(27r|x|) _1 and J p = A. We have the following bounds on the potential p*\. 



(P*\.n(*)> I ttttit^- 



(i) (Lower bound) For all x G 1R 2 \{0} ; 

max{|x|, \y\} 

(ii) (Upper bound) For all x G M 2 \{0}, 

(p*\. (a) < 2V2\\x\- 1/2 + 3. 

Moreover, for any 5 > there exists R = R(p, 5) > and a universal constant C\ > 
such that for any \x\ > R, 

/ i i i \ / \ ^ A ~\~ $ „ ln(b|) f 

(p * I • I )(%) < + Ci—^ 1 / p(y)dy. 



3\x\/2>\y\>\x\/2 
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3 Thomas- Fermi Theory 



In this section, we introduce the two-dimensional Thomas-Fermi (TF) theory which will 
turn out to be the main tool to understand the ground state energy and ground states. The 
three-dimensional TF theory was studied in great mathematical detail by Lieb-Simon [T5|l9]. 
In fact, the simplest version of TF theory (see [T2], Chap. 11) is sufficient for our discussion 
here. 

Definition 13 (Thomas- Fermi functional). For any nonnegative function p G L X (IR 2 ) we 
define the TF functional as 

S TF (p) := J (^ p \ x )-E^ + (An)- 1 [\x\- 1 -l}l^dx + D(p). 

R 2 

For any A > we define the TF energy as 

E TF (X) := inf {s TF (p)\p > 0, ||p|| xl(R2) < A} . (5) 

Remark. (i) The term irp 2 comes from the semiclassics of the kinetic energy while — J p(x) \x\ 
and the direct term D(p) = ~ ff p(x)p(y)\x — y\~ 1 dxdy stand for the Coulomb inter- 
actions. 

(ii) The appearance of (47r)~ 1 [|x|~ 1 — 1]+ ensures that the TF functional is bounded from 
below. In fact, 

£ TF (p) = J -[P^)~^)dx+ J (« p \ x )- P -^d x + D{p)-\ 

\x\<l ' \x\>l 



> -Ip-- 4 - 

(iii) If pJ F is the ground state of the above TF theory then Zp^ F is expected to approximate 
the density p^, NZ of a ground state ^n,z of H N Z with iV ps XZ (in some appropriate 
sense). In other words, the Z-dependent TF theory can be defined from the above TF 
theory by the scaling p i— > Zp. 

Basic information about the TF theory is collected in the following theorem. 

Theorem 14 (Thomas- Fermi theory). Let A > 0. 

(?) (Existence) The variational problem (j5J) has a unique minimizer p^ F . Moreover, the 
functional A i— >■ E TF (\) is strictly convex, decreasing on (0, 1] and E TF (X) = E TF (1) if 
A > 1. 

(ii) (TF equation) p^ F satisfies the TF equation 

2npl F (x)=[\xr-(pl F *\.r)(x)-pJ F ] + 
with some constant pJ F > if A < 1 and pJ F = if A > 1 . 



S 



(iii) (TF minimizer) p^ F is radially symmetric; j p^ F = min{A, 1} and 

< \x\-~ 1 - 2npJ F < C\x\- l/2 for all x ^ 0. 
Moreover, supp pJ F is compact if and only if A < 1 . 

Remark. Henceforth we shall always denote by C some finite positive constant depending 
only on A > (the total mass in the TF theory). Two C's in the same line may refer to two 
different constants. 

Proof, (i-ii) Formula implies that p i— > S TF (p) is strictly convex. Therefore, the existence 
and uniqueness of the TF minimizer, and the TF equation follow from standard variational 
methods similarly to the three-dimensional TF theory (see [T2], Theorems 11.12 and 11.13). 
The property of pJ F is a consequence of the TF equation and is shown in Lemma [15] below. 

That E TF (X) is decreasing follows from the definition. When A > 1, E TF (X) = E TF (1) 
since pJ F = pJ F (by (iii)). When A G (0, 1], the TF energy is also strict convex because the 
unique minimizer satisfies J pJ F = A (by (iii)) and the TF functional is strict convex. 

(iii) Since the TF functional is rotation invariant and the minimizer is unique, it must 
be radially symmetric. The inequality < — 27rpJ F < C|x| -1 / 2 follows from the TF 
equation and the following estimate in Lemma [12], 

(p TF *|.r 1 )(x)<2V2A|x|- 1 / 2 + 3. 

We defer the proof that f pJ F = min{A, 1} and property of supp pJ F to Lemma [TBI □ 



3.1 Thomas- Fermi Equation 

Lemma 15 (TF equation). Assume that p is a nonnegative, radially symmetric, integrable 
solution to the TF equation 

2K P (x)=l\x\- l -(p*\.\- 1 ){x)-p] + (6) 

for some constant p > 0. 

(i) If p > then j p < 1 and supp p is compact, 
(ii) If p = then J p = 1 and 



/ p(x)dx > e 2 ^ for all r > 0. 

J\x\>r 



Proof. Denote j p =: A > 0. For r > we shall write p(r) instead of p(x)\\ x \ =r . 

1. We start by proving A < 1. Since p is nonnegative and radially symmetric, we have 
by Lemma fT2l 

p(y) 



(p*I-I _1 )^)> 

Hence, the TF equation ([6]) yields 

27rp(a;)|x| < 





x\p(y) 


max{ \x 


•> 


y\} 



R2 max{|x|, 1 



dy — p\x\ 



dy. 



for all s^0. 



(7) 
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For any e G (0, A), we can find R £ > such that f\ x \ >R p = £■ When \x\ > R £ , using 

V max{|x|, \y\} J\ y \<R e 

we can deduce from ([7j) that 

2ttp(x)\x\ < [1 — A + e — p|x|]+ < [1 — A + e — pR £ } + for all \x\ > R £ . 

Since f\ x \ >R p = £ > 0, there exists \x\ > R £ such that p(x) > 0. Therefore, it follows from 
the latter estimate that 

1 - A + e - pR £ > for all e G (0, A). (8) 

For any p > 0, (JSJ) implies that A < 1. 

2. If /i > then (J8J yields 

limsupi? e < i?o : = P (1 — A). 

Since Jj^^ P = £ and lim sup e _> i? £ < i?o, we get JL^j^P = 0. Thus suppp C {\x\ < Ro} 
and A < 1 (because R > 0). 

3. From now on we assume that p = 0. We shall prove that A = 1. Suppose that 
A < 1 — 3e for some e > 0. Because p is nonnegative, radially symmetric and p(x) < (27r|a;|) _1 
(due to the TF equation ©), by Lemma [12] we can find R > and Ci > such that 



(p * I ■ T 1 )^) < \^ + Ci^^ / p{y)dy for all |x| > R. (9) 

pi pi «/3|a:|/2>|i/|>|3;|/2 

Define E\ :— e/C\ and 



A := \ r > R : / p(y)dy < ' ' 



'3r/2>\y\>r/2 l n ( r , 

If \x\ G A, then ([9]) gives {p * \ ■ \^ 1 )(x) < (1 — and the TF equation fl6]) with p = 

gives 



2np(x) 



x 



(p * |s| x )(x) 



> — 

a; 



Taking the integral of the previous inequality over {x G M 2 : |x| G A} one has 
oo > 2tt j p > 2tt y p(x)rfs > y ^rfx = 2ixeC 1 (A) 



\x\&A \x\eA 



where C l is the one- dimensional Lebesgue measure. Thus A has finite measure, and conse- 
quently we can choose a sequence {Rn}™^ C M\A such that 3R < 3R n < R n +i < 4i? n for 
all n > 1. Because R n > R and i? n ^ A we have, by the definition of A, 



/ 



3i?„/2>|y|> J R„/2 
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p(y)dy > ; — ■— for all n > 1. 
vy; ln(i? n ) 



Taking the sum over all n G N and using R n+ \ > 3R n , we find that 



p > Yl / p(y) d v > Yl 

n=l„r, ,„^T r, ,„ n=l 



oo 



ln(K) 

± 3i?„/2>|y|>iJ n /2 " " 

On the other hand, since i? n+ i < AR n we get < A n R l < [4(1 + #i)] n for all n > 1. 
Therefore, 

oo OO ^_ 

rnTi?J " ^ nln(4(l + iZj) = + °° 

The last two inequalities yield a contradiction. 

4. Finally, we show the lower bound on J^ >r P- With fi = and A = 1, inequality (J2J) 
becomes 

2irp(x)\x\ < J (l - j^jj p{y)dy for all x ^ 0. (10) 

l«l>M 

Denote 



0(7 J 



/ ( x ~ Si ) = 27r / ^ s ~ ^ p ^ ds - 

\y\>r r 

Then g(0) = 1, g(+oo) = and 

oo 

g'(r) = —2n J p(s)ds < 0, g"(r) = 2ixp{r) for all r > 0. 



Thus ( flop can be rewritten as 

rg"{r) < g{r) for all r > 0. 
Note that go(r) := e -2 ^ satisfies go(0) = 1, (?o(+oo) = and 

r<(r)-^ (r) = ^=e- 2 ^>0. 

Therefore, h(x) := gr(a;) — ^o(^) satisfies that h(0) = h(+oo) = and rh"(r) < h(r). If the 
set U := {r > : /i(r) < 0} is not empty, then h is a strict concave function on this open 
set. By the maximum principle and h(0) = h(+oo) = 0, we can argue to get a contradiction. 
Thus h(r) > for all r > 0. This yields f,, >r p(x)dx > g(r) > go(r) = e~ 2 ^. □ 

4 Semiclassics for the TF Potential 

In this section, we consider the semiclassics for the TF potential 

n TF (x):= \x\- 1 -(pl F *\.\- 1 )(x)-fiJ F . 
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From the TF equation and the properties of the TF minimizer (see Theorem IT4"|) we have 
[V A TF ]+ G L 1 ^ 2 ) and 

|^ a tf 0e) - la;]- 1 ! < C(\x\- 1/2 + 1) for all x ^ 0. 

The following theorem will turn out to be the main ingredient to prove Theorems [T] and 
El The parameter h will eventually be replaced by {2Z)~ l l 2 in our application. 

Theorem 16 (Semiclassics for TF potential). When h — > + one has 

Tr[-h*A-V?*]_ = -(S^W ([V^(x)} 2 + - [\x\^ - I] 2 ,) dx 

Jr 2 

+ (4h 2 )- 1 [\n(2h 2 ) + c H ] + o(h- 2 ). (11) 

where c H = -3 ln(2) - 2-f E + 1 « -2.2339. 

Moreover, there is a density matrix 7^ such that 

Tr [{-h 2 A - V? F h h ] = Tr [—h 2 A - V? F ] _ + o(^ 2 ) (12) 

and 

2h 2 Tr( 7h ) < | pJ F , D((2h 2 )p lh - pD = o(l). (13) 

Note that fill I) is a special case of Theorem [3J In this section, we shall prove fllip in 
detail. The proof of Theorem [3] is provided in the next section. 

As in [26] we shall prove the semiclassical approximation (fTTj) by comparing with the 
hydrogen. In fact, because of the hydrogen semiclassics ([2]), the approximation fTTT]) is 
equivalent to 

Tr [-h 2 A - V? F ] _ - Tr [-h 2 A - \x\- x + l] _ 
= -(Snh 2 )- 1 [ {[Vni-ll^-^Ddx + oih- 2 ). (14) 

</R 2 

4.1 Localization 

To treat the singularity of the TF potential we shall distinguish between three regions. In 
the interior region (close to the origin), we shall compare directly with hydrogen; while 
in the exterior region (not too close and not too far from the origin) we can employ the 
usual semiclassical techniques; and finally, the region very far from the origin has negligible 
contribution. 

Definition 17 (Partition of unity). Let ip be a nonnegative, smooth function (with bounded 
derivatives) such that ip(x) = 1 if \x\ < 1 and (p(x) = if \x\ > 2. Choose r := h 1 ^ 2 , 
A := | In h\ and denote 

$i(x) = <p(x/r), 

$ 2 (x) = (l- V 2 (x/r))^ 2 ^x/A), 

$ 3 (x) = (l-^Or/A)) 1 / 2 . 

Then YLi ^ = 

supp$i C < 2r}, supp $2 C {r < |x| < 2A}, supp$3 C {|a;| > A}. 
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The localization cost is controlled by the following lemma. 

Lemma 18 (Localization). Let V be either V^ TF or — 1). When A = \ lnh\ and 

r = h 1 ! 2 — > + one has 

Tr[-/i 2 A-V]_ = ^Tr[$ i (-/i 2 A-y)$ i ]_ + o(/i- 2 ) 

i=l,2 

Note that in the sum on the right-hand side the contribution of region supp $3 does not 
appear. 

Proof. 1. To prove the lower bound, using the IMS formula 

3 3 
-A = ^(" A " u )®i with u ■= I V$i ' 2 - Cr_2l {l-I<2A} 

i=l i=l 

one has 

3 

Tr[-/i 2 A - V]_ > ^Tr[$ ?; (-/i 2 A - V - Ch 2 r- 2 l {[x[ < m )$i]-. 
i=i 

The term involving $ 3 has negligible contribution. Indeed, since supp $3 C {\x\ > A}, it 
follows from the Lieb-Thirring inequality ([3]) that 

Tr[$ 3 (-/i 2 A - V - C/i 2 r- 2 l { | iC |< 2A} )$ 3 ]- 

> Tr[-h 2 A - 1 { |,|>A } (V + + C7rV- 2 l {N < 2A} )]_ 

> -L 1>2 h~ 2 / [V + {x) + Ch 2 r- 2 l {lxl<2A} ] 2 dx = o(h- 2 ). 

J\x\>A 

Here note that limA-s.00 f\ x \>& V 2 = since 1{>|>i}V + G L 2 (IR 2 ) and /i 4 r~ 4 Jj x | <2A — > because 
A=\\nh\. 

Moreover, for % = 1, 2, if we denote 

7, := l(_oo,o] {M~h 2 A -V- Ch 2 r- 2 )§ t ) 

then Tr($i7i$i) < Ch~ 2 {\ \ah\ + A 2 ) by Lemma [8] (ii) . Therefore, 

Tr[$;(-/i 2 A- V -Ch 2 r- 2 )<5>i]_ = Tr[$;(-/i 2 A - V - Ch 2 ^ 2 )®^] 

= Tr[^(-/i 2 A - V)$ l7l ] - Ch 2 r- 2 Tr($ i7 i$i) 
> Tr[$,(-/i 2 A- V)$,]_ +o(/T 2 ). 

2. To show the upper bound, we choose 

7 W := l(-oo,o] (H~h 2 A - V)$0 , 7 (0) := ]T ^7<$i- 

i=l,2 

Since < 7® < 1 (i = 1, 2) and E;=i, 2 < 1 we have < 7 (0) < L Thus ' 
Tr[-/i 2 A - V]_ < Tr[(-/i 2 A - V) 7 (0) ] = ^ Tr[$ ?; (-/i 2 A - 

i=l,2 

□ 
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4.2 Hydrogen Comparison in Interior Region 

In the interior region, we shall compare the semiclassics of the TF potential directly with 
hydrogen. Note that 



< CrbT 2 = o{h~ 



because \V^ F — \x\ 1 \ < C(\x\ l l 2 + 1) and supp $1 C {\x\ < 2r}. This inequality is the 
semiclassial version of the following bound. 

Lemma 19 (Hydrogen comparison in interior region). When r = h 1 / 2 — > we have 

Tr [$! (-h 2 A - V? F ) - Tr [$ x (-h 2 A - {x^ 1 + l) $J_ = o(h~ 2 ). 

Proof. The lower and upper bounds can be proved in the same way. We prove for example 
the upper bound. If we denote 

7 {1) := l(-oo,o] [$i {~h 2 A - Ixl' 1 + 1) $J 

then by Lemma [8] (ii), 

Tr[$i 7 (1) $i] < Crh~ 2 \ lnh\ 1/2 . (15) 
By using |V A TF (:r) — |a;| _1 + 1| < C(\x\~ 1/2 + 1) < Cr~ 1/2 for x e supp $i we get 

Tr [$! (-/i 2 A - Is]" 1 + 1) $i] = Tr [$ x (-/i 2 A - |ac| -1 + l) $i 7 {1) ] 

> Tr [$! (-/i 2 A - V A TF ) $i7i] - C r_1/2 l^iT^i] 

> 'IV[$ 1 (-/ l 2 A-'^ F )$ 1 ]_ + o(/i- 2 ). 

□ 

4.3 Semiclassics in Exterior Region 

In the exterior region, the standard semiclassiccal technique of using coherent states [HJ EI] 
(see also 112, 241) is available. 



Definition 20 (Coherent states). Let g be a radially symmetric, smooth function such 
that < g(x) < 1, g(x) = if |x| > 1 and f R2 g 2 (x)dx = 1. For s > (small), denote 
g s (x) = s^ 1 g(x/s) and 

n s , n ,p = |/ a> uj») where f s , u , P ( x ) = eiP ' X 9s{x ~ u) for all u,p E R 2 . 

From the coherent identity, 

(2tt)" 2 J I U S)UiP dpdu = I on L 2 (R 2 ), (16) 

l 2 xl 2 
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it is straightforward to see that for any density matrix 7 and for any potential V satisfying 
V+ E L^R 2 ), 

Tr [-/i 2 A 7 ] = (2n)- 2 JJ Tr [-/i 2 An s „ 7 ] dpdu 

= (2tt)- 2 J J /iVTr[n SitiiP 7]^-||V^||i 2 /iV 2 Tr(7), (17) 

Tr[(-V * g 2 s )i] = (2tt)- 2 / j Tr[(-V * g 2 )U s ^ pl ]dpdu 

= (2tt)- 2 J J -V(u)Tr[H StU>pl ]dpdu. (18) 

Motivated by ffTB"]) . it is useful to have some estimate for (V — V * g 2 s ). The proof of the 
following lemma can be found in the Appendix. 

Lemma 21. IfV is either V^ F or (jx)" 1 — 1) and A = \ \nh\, r = h 1 ! 2 , s = h 2 ^ 3 then 

\V - V * g]\ 2 (x)dx <Ch l '\ 



r<\x\<2K 

Lemma 22 (Semiclassics in exterior region). Let V be either V^ F or (|x| _1 — 1). When 
A = I lnh\ and r = h 1 / 2 — > one has 

Tr [$ 2 (-h 2 A - V) $ 2 ]_ = -{Znh 2 )- 1 J V}(x)$l(x)dx + o(/T 2 ). 

Proof. 1. To prove the lower bound, we choose the density matrix 

72 := l(-oo,o] [$2 {~h 2 A - V) $ 2 ] . 

Taking s = h 2 / 3 and using identities f fT7|) and (TT~8]) we can write 

Tr [$ 2 (-/i 2 A-y) $ 2 ]_ = Tr [(-h 2 A - V) $ 272 $ 2 ] 

= (2tt)- 2 ^ [/iV-y(n)]Tr[n s , UiP $ 272 $ 2 ]^ 

+ Tr [(V * £ 2 - V - Ch 2 s' 2 ) $ 272 $ 2 ] . (19) 

2. To bound the second term of the right-hand side of (119)) . we can apply Holder's 
inequality, Lemma [8] (i) with f2 := supp$ 2 C {r < |a;| < 2A} and Lemma I2T1 to get 



> -\\ v * 91- v \\t.voAp*™*2\\l*<w\ -Ch 2 s 2 Tr[$ 272 $ 2 



Tr [(V * g 2 - V - CTiV 2 ) $ 272 $ 2 ] 

> -Ch- 2 \\V * g 2 - V\\ L2{Q) \\V + \\ L2(n) - Cs- 2 \\V + \\ L2{n) M" 2 

> -Ch^h^hih^ 2 -Cs- 2 \\nh\ l ' 2 K = o(h- 2 ). (20) 



For the first term of the right-hand side of (fT9l . because 

< Tr [n SiUjP $ 272 $ 2 ] < Tr [IT SiUiP $ 2 ] = ($ 2 * g 2 )(u) 
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we obtain 



(2tt)- 2 J I [h 2 p 2 - V(u)] Tr [n s , UiP $ 272 <l> 2 ] dpdu 
> -W 2 jj [hV-V(u)]_(^l*g 2 s )(u)dpdu 
= -(Svr/^r 1 ! V 2 (u)($ 2 *g 2 )(u)du 

= -(87th 2 )- 1 Jv 2 (u)$i(u)du + o(h- 2 ). (21) 



Here the last estimate follows from 

I Kj(u)|$| - $2 * £ s 2 |(w)dw < Csr- 1 [ V + {u) 2 du < Csr- 1 ] lnr| = o(/T 2 ), (22) 

J J\x\>r/2 

where we have used |$ 2 — ^2 * 5s I ( x ) — C sr_1 l{M>r/2} when |x| > r ^> s. Replacing ( 12"U|) 
and ()2ip into ffT9l) we get the lower bound in the lemma. 
3. To show the upper bound, we choose 

7® := (2tt)- 2 J J U SjU>p dpdu, M := {(u,p) : /^V - < 0} . 

M 

Using the coherent identity ({TBI) and the IMS formula, it is straightforward to compute that 
Tr [$ 2 (-h 2 A - V) $ 2 ]_ < Tr [$ 2 (-/i 2 A - V) $ 2 7 (2) ] 
= (2tt)- 2 / / (/ S)UjP (x)| $ 2 (z) [-h 2 A x - V(x)) <f> 2 {x) \fs,u, P { x )} L z(WL2,te) d P du 



M 

,h 2 . /i 2 



(2vr)- 2 / / (e^| - (^ s $ 2 ) 2 yA - y A(^$ 2 ) 2 + h 2 \V(g s $ 2 )\ 2 - (g s $ 2 ) 2 V \e**) dpdu 



M 



(2n)- 2 // hY(^2*9t)(u) + h' / \V(g s ^ 2 )(x)\ 2 dx-((^ 2 V)*gi)(u) dpdu 



M 



= -{^h 2 )- 1 J V 2 (u)$ 2 2 (u)du + (Snh 2 )- 1 J [Vl^l * g 2 ) - V 2 ® 2 ] du 

+(4tt)- 1 jj V + (u)\V(g s $ 2 )(x)\ 2 dxdu+ (Anh 2 )- 1 J $ 2 U [V+ - (V+ * g 2 )] du. (23) 

4. Finally we verify that the last three terms of the right-hand side of (1231) are of o(h~ 2 ). 
The second term was already treated by f l22l) . Using 



/ \V(g s $ 2 )(x)\ 2 dx<C(r- 2 + s- 2 )l { 

r/2<|«|<3A}- 

JR 2 



we can bound the third term as 



V + (u)\V(g s $ 2 )(x)\ 2 dxdu < C(r- 2 + S - 2 ) / V + {u)du = o(/T 2 ). 

Jr/2<\x\<3A 
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To estimate the last term, we introduce a universal constant A > such that V(x) > 
when |x| < 2A (such A exists since \V(x) — < C(|x| _1 / 2 + 1)). Using V + (u) = V(u) 

and (V + * g 2 s ){u) — (V * g 2 s ){u) when \u\ < A , and Lemma I2T1 we get 

/ \$lV\.\V + - V + * gl\du < \\V\\ L2{n) \\V * g 2 s - V\\ < C\ \nh\ l ^ s = o(l) 

J\u\<A V ; 

where Q = supp $ 2 C {r < \u\ < A}. On the other hand, because |U(u)| < C when \u\ > A 
and V + G L^R 2 ), 

/ \$lV\.\V + - V + * gl\du < C \\V + * g\ - V+ \\ L 1(R2) = o(l). 

Thus the last term of the right-hand side of f[23l) is also of o(h~ 2 ). This completes the 
proof. □ 

Lemmas [TBI dH and [25] together yield ( TT4"]) , which is equivalent to (ITT]) . 
4.4 Trial Density Matrix 

The last step in proving Theorem [16] is to construct a trial density matrix. 
Lemma 23. There exists a density matrix 7^ satisfying ([12]) and (TIB"]) . 



Proof. Recall that we always choose A = | lah\, r = /i 1//2 and s = h 2 ^ 3 . 

1. From the proof of Lemmas [18] [19] and [22], if we choose the density matrices 



« := 1(^,0] [$! (-/i 2 A - l^r 1 + 1) $J 



7 



7 (2) 



(2tt) 2 ^ Il s ^ p dpdu, 



h 2 p 2 -V^ F (u)<0 



7 (0) . = $ l7 (l)$ 1 + $ 27 ( 2 )$ 2 



then 



Tr[-/i 2 A - U A TF ]_ = Tr[ (_/,2 A _ vf F ) 7 (°)] + o(/i- 2 ). (24) 



2. Using the coherent identity f[16]) and the TF equation p^ F = (2ix) 1 [U !v TF ] + , we can 
compute explicitly that 



P^){x) := ^ 2 \x,x) = (2tt)- 2 



yy n SiUiP (x, x) dpdu 

V? F (u)<0 

{^r\[vn + * £){x) = (2h 2 r\ P r * ^ 



7 

/l 2 p 2_yTF( u )<Q 

2\-l/rT/TFl , „2\/ \ /oi,2\-l/„ TF , „2 



Therefore, 

2/i 2 p 7(0 ) = 2/i 2 p $i7(1)$i + $ 2 pJ F * <? 2 . (25) 
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Since / pJ F * g 2 s = J pJ F and Tr^W^] < Crh- 2 \ \nh\^ 2 (see (USD), we have 

2h 2 [ p^(i )01 {x)dx< [ pl Y {x)dx + Cr\\nh\ 1/2 . (26) 
Jr 2 Jr 2 

On the other hand, we can write from (I25p that 

2/i 2 p 7 (o) - pI F = 2/i 2 p* l7 w* 1 + $ 2(pI f * 9 2 s - p7) + (i - $2)pI f - 

Since pJ F e L 4/3 (M 2 ), we have pJ F * g 2 - pJ F and ($| - l)p^ F converge to in L 4 / 3 (M 2 ). 
Moreover, using Lemma M we have 2/i 2 p $l7 (i) $1 — > in L 4 / 3 (IR 2 ). Thus 2h 2 p 1 ( ) — pJ F — > 
in L 4 / 3 (R 2 ). Since the Coulomb norm is dominated by the L 4 ' 3 -norm (see Theorem [TT]) , 
we then also have 

D(2h 2 p y{0) - p T x F ) -> 0. (27) 

3. Finally, we choose £ such that |ln/i| _1 ^> £ r | In /z.1 1 / 2 (e.g. £ = r 1 / 2 = h l / A ) and 
define 

lh := (l-£) 7 (0) . 
Then using (12"6"|) and £ ^> r| In /i| 4 / 2 we have 



2/i 2 Tr( 7ft )<(l-£)(l + Cr|ln/i| 1 / 2 ) / P^ < / 



pI f 



for h small enough. Moreover, since \lnh\ 1 ^> £, the inequalities f[24"|) and (I27p still hold 
true with 7^ replaced by 7^. □ 



5 Proofs of the Main Theorems 
5.1 Ground State Energy 

Having the semiclassics in Theorem [161 the proof of Theorem [1] is standard (see [9]). 

Proof of TheoremUi 1. We first prove the lower bound. Taking any (normalized) wave 
function \& G /\ i=1 L 2 {R 2 ), we need to show that 

#iv,^) > In Z + £ TF (A)Z 2 + o(Z 2 ). 

Starting with the Lieb-Oxford inequality (Theorem [9]) 

E F^^l >^(P*)-^LO / pf , 
V l<i<j<N / ^ 

we want to bound / p^ 2 . It of course suffices to assume that (^>,Hn^^) < 0. Using the 
Lieb-Thirring inequality fl4]), the hydrogen spectrum in Theorem [5] and Tr(7*) = N < CZ 
we arrive at 

> 4 (*, H NjZ ^) > Tr[-A 7 *] + Tr [(-A -AZ^ 1 )^] 

> K 2 [ p%{x)dx-CZ 2 \\nZ\ 
Jr 2 
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By Holder's inequality and J pq, = N < C Z again we conclude 



Py(x)dx < I I p\{x)dx 



Thus the Lieb-Oxford inequality gives 



1/2 



r \ 1/2 

/ p^{x)dx) <CZ^ 2 \\nZ\ 1 ' 2 . 



(^,H N ^) > Tr 



—A-Z\x 
2 1 



-i 



7* 



+ D(p^)-CZV 2 \lnZ\ 1 



/2 



ZTr [(-^A - V? F ) 7w ] - Z 2 [p7{N/Z) + D{p 
+Z 2 D{Z- l p^ - pl ¥ ) - CZ 3 / 2 \ InZ] 1 / 2 . 



TF 
A 



(28) 



For the lower bound, we can ignore the nonnegative term D(Z 1 p l9 — p^ F ) — 0- With the 
semiclassics of the TF potential in Theorem [TBI and h 2 = (2Z)^ 1 , one has 



Tr [(-{2Z)-^ - VD 7*] > Tr [-pZ^A - V? F }. 



> — Z\nZ 

~ 2 



-(4-r 1 / {[vr{x)t-[\xr-i] 2 + )dx + y 



+ o(Z). 



Together with N/Z — > A, we obtain from (128]) that 



where 



(tf , tf^z*) > ~\z 2 In Z + e(A)Z 2 + o(Z 2 ) 



<A) := -(4VT)- 1 1 ([^(x)] 2 , - [Ixr 1 - l] 2 + ) dx - pJ F X - D(pD + 1 C H 



By the TF equation 2irpJ = [V A ] + we have 



^2r /TFi2 



^HN^-pr*!-!- 1 -^]- 



Replacing this identity and pJ F X = pJ F J pJ F into the definition of e(A), we see that e(A) = 
E TF (X) + c H /2. Thus we get the lower bound on the ground state energy. 

2. To show the upper bound, because A h-> E tf (X) is continuous, it suffices to show that 
for any < A' < A fixed, one has 

E(N, Z) <-^Z 2 \nZ + E TF (X')Z 2 + o(Z 2 ). 

Using Lieb's variational principle (see Theorem [TUl) we want to find a density matrix 7 such 
that Tr(7) < N and 



Tr 



-A - Zld -1 
2 1 1 



7 



+ D( Pl ) <--Z 2 \nZ + E Lt (X')Z 2 + o{Z A 



This condition can be rewritten, using the same calculation of proving the lower bound (see 
(|28|)), as 

Tr [(-{2Zy 1 ^ ~ V7) 7] + ZD(Z' l p, - p T x F) < Tr [-(22)^ - V? F }_ + o(Z). (29) 

According to Theorem [TBI with h 2 = (2Z)~ 1 , we can find a trial density matrix 7 satisfying 
such that Tr( 7 ) <Zj pJF < X'Z. Since N/Z A > A', one has Tr( 7 ) < X'Z < N for 
Z large enough and it ends the proof. □ 
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5.2 Extensivity of Neutral Atoms 

Proof of Theorem^ Let 9r be a smooth function such that 9r(x) — if \x\ < R and 
9r{x) = 1 if \x\ > 2R. From the proof of Theorem [16] and Theorem [TJ, we have, with 
7 := 7*iv, z and /i 2 = (2Z)' 1 , 

Tv[(-h 2 A - ^ TF ) 7 ] = Tr[-/i 2 A - ^ TF ]_ + o(h~ 2 ). 

Using the localization as in Lemma [18] and the semiclassics of Lemma [22] we get 

Tr[9 R (-h 2 A- V? F )9 Rl \ < Tr[0 jR (-/i 2 A - V^ F )9 R ]- + o{h~ 2 ) 

= -(871k 2 )- 1 [ [V? F ] 2 + (x)9 2 R (x)dx + o(h~ 2 ). (30) 



On the other hand, since V± < \x\ 1 < R 1 in supp6^, 

Tr[^(-/i 2 A - V? F )9 Rl ] > -iT 1 Tr[^ T ^] = -iT 1 / 9 2 R (x)p^(x)dx. 
Putting (130]) and (131]) together we arrive at 

^ 2 (x)p 7 (x)dx > R{8nh 2 )- 1 J [V? F \\(x)9 2 R (x)dx + o(/i- 2 ). 
Replacing h 2 = (2Z) _1 , we can conclude that 

p^(x)dx > / 9 2 R (x)p 1 {x)dx > CrZ + o(Z) 



\x\>R 



where 



C R := ii!(47r)- 1 ( [V? F ] 2 + (x)9 2 R (x)dx > ttR [ {p TF {x)) 2 dx. 

J J\x\>2R 



(31) 



Note that Cr > because supp pJ F is unbounded (see Theorem [TJJ . □ 



5.3 Semiclassics for Coulomb Singular Potentials 

Proof of Theorem 0. We shall show how to adapt the proof of (ITT]) in the previous section 
to the general case. We however leave some details to the readers. By scaling we can assume 
k = 1. 

1. The main difficulty of the general case is that we do not have the estimate in Lemma [211 
in the exterior region. Therefore, we need a more complicated localization. Let r = h 1 ^ 2 , s = 
h 2 / 3 and let g s be as in Definition [20] For any e > small, denote 

W{e, h) := [ \V\ 2 dx [ \V+ - V + * g 2 \ 2 dx 
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Because V E L 1 2 OC (R 2 \{0}), for any s > fixed we have W(e, h) — >■ as h — >■ + . Therefore, 
we can choose e = e(h) such that e(/i) > | \nh\~ 1 , e(h) — > and W(e(/i), /i) — > as /i — > + . 
Let <y? as in Definition [T7] and define 

$i(x) = (f{x/r), 

$ 2 (x) = (l-<p 2 (x/r))VMx/e), 

$ 3 (x) = (l-^x/ef'Mxeffl, 

$ 4 (x) = (1 - ^(xe/2)) 1 / 2 . 

Then £)f =1 = 1, supp$i C < 2r}, supp $ 2 C {r < |x| < 2e}, supp $ 3 C {e < \x\ < 
e' 1 }, and supp $ 4 C {\x\ > (2e)~ 1 }. 

2. Following the proof of Lemma [TBI we can show that 



Ti[-h 2 A -V]-=J2 Tr[$i(-/i 2 A - V)$i]- + o(/T 2 ) 



(32) 



i=l 



Note that the assumptions l-oi>i}^+ £ L 2 (IR 2 ) and \ \nh\ > e^ 1 — > oo is sufficient to bound 
the contribution of the region supp$ 4 by the Lieb-Thirring inequality (|2])- To control the 
localization cost in the region supp $3, we may use Lemma [8] (i) instead of Lemma [8] (ii). 

3. Because \V(x) — | a; | 1 + 1| < C(|a:|~ + 1) < Cr~ e for x G supp$i, we can follow the 
proof of Lemma [19] to get 



Tr 



$1 (-h 2 A - V) $1 



Tr 



$1 (-h 2 A 



x 



-1 



+ l)$i +o(h' 



(33) 



4. Adapting the coherent state approach in the proof of Lemma [221 we can show that 



Tr 



$ 3 (-h 2 A - V) $ s 



2\-l 



V 2 (x)$ 2 3 (x)dx + o(h- 2 ). 



(34) 



To obtain the lower bound it suffices to consider Tr[$3 (—h 2 A — V + ) $3]- and then use the 
assumption W(e(h), h) — > instead of Lemma I2T1 in (120]) . When proving the upper bound, 
the assumption W(e(h), h) is again enough to estimate the last term of f[23|) . 
5. In the intermediate region supp $2 C {r < \x\ < 2e}, we have 

V 1 (x) : = (a;]" 1 + C\x\~ e > V{x) > \x\-~ 1 - C\x\~ e =: V 2 (x) > 0. 

We start with the lower bound 



Tr 



$2 (-h 2 A - V) $ 2 



> Tr 



$2 (~h 2 A - $ 2 



Using the coherent state approach as in the proof of Lemma [221 we can show that 



Tr 



V! 2 (x)$2(x) 2 rfx + o(/i- 2 ). 



$ 2 (-h 2 A - Vr) $ 2 _ = 
To do that, we just need to replace Lemma I2T1 by the following estimate 

\Vx-Vx * g 2 \ 2 dx < Cs 2 |lnr|. 

r<\x\<2e 
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Moreover, since supp $ 2 C {r < \x\ < 2e} with e = e{h) — > and \V — Vi\ < C\x\ e , we 
have 

r V?(x)$ 2 (x) 2 dx = J V 2 {x)§ 2 {x) 2 dx + o{fr 2 ). 
Therefore, we arrive at 



$ 2 (-h 2 A - V) $ 2 > / V 2 (x)$ 2 (x) 2 dx + o(h 



Tr 



Similarly, again using the coherent state approach we get the reverse inequality 



Tr 



$ 2 (~h 2 A - V) $ 2 < Tr 

V 2 2 {x)$ 2 {x) 2 dx + o{h- 2 ) = 



$ 2 (-h 2 A - V 2 ) $ 2 ^ 
V 2 (x)$ 2 (x) 2 dx + o(h~ 2 ) 



Thus we obtain the semiclassics 



Tr 



$2 (-h 2 A - V) $ 2 



V 2 (x)$ 2 {x) 2 dx + o(h~ 



6. The desired semiclassics follows from (|32l . (j33J), fl34|) and ([35 



(35) 
□ 



A Appendix 

In this appendix we provide several technical proofs. 

Proof of Theorem^ (i) The HVZ Theorem indeed holds for all dimension d > 2 (see e.g. 
[8] Theorem 2.1 for a short proof). The decay property is essentially taken from [18] where 
the only change is of solving equation (3.8) in [18]. In fact, the two-dimensional solution 
w 2 (r) (with r = \x\) is obtained by scaling the three-dimensional solution w^r) in [IB] as 

W 2 = Ws\ £ ^4 :£t Z^4:Z,r^r/2. 

(ii) The proof of Zhislin's Theorem is standard and there is no difference between two 
and three dimensions. The idea is that by induction we can use the ground state Hn,z to 
construct a (N + l)-particle wave function with strictly lower energy whereas N < Z. It 
should be mentioned that some certain decay of the ground state is necessary to control the 
localization error when we consider the cut-off wave function in a compact set. 

(iii) The asymptotic neutrality follows from the original proof in three dimensions of Lieb, 
Sigal, Simon and Thirring [T3]. The key point of their proof is the construction of a partition 
of unity. But a partition of unity in three dimensions obviously yields a partition of unity in 
two dimensions, hence this part of the proof can be adopted. Note that the Pauli exclusion 
principle enters when solving the hydrogen atom. □ 



Proof of Lemma For any m G N, one has 
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Using Euler's approximation 

m | 

^ V+lj2 = + 2 lll(2) +lE + 

we get 
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which implies ([Q). Moreover, fl2]) follows from ([[]) by scaling x h-> (2/i 2 ) 1 x, namely 

1 



Tr [-/i 2 A - Ixl" 1 + //]_ = (2/i 2 )- 1 Tr 



-A - Ixl" 1 + 2/i 2 /i 



Proof of Lemma\^ (i) For any constant a > 0, using the Lieb-Thirring inequality ( 
have 

> Tr[(-/i 2 A - V)0 7 0] > Tr[(-(/i 2 /2)A + a)0 7 0] + Tr[(-(/i 2 /2)A - (a + ^+).l n 



□ 

we 



> Tr[(-(/i 2 /2)A + a)0 7 0] - AL h2 h~ A \\a + V + \\\ 2 
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Choosing a = 1 and using Tr[— A070] > and (1 + V + ) G L 2 {VL) we get Tr[070] < 00, 
namely 070 is trace class. On the other hand, choosing a = and using the Lieb-Thirring 
inequality (jlj) to estimate Tr[— A070], we arrive at 



ffa(x)dx < Chr*\\V + \\ 



Because supp p^ C Q, the above estimate and Holder's inequality yield the desired bound 

011 J pft/cp for an y a e [M]- 

(ii) We can use the same idea of the above proof. The only adaption we need in this case 
is to use both of the Lieb-Thirring inequality ([3]) and the hydrogen semiclassics fl2]) to bound 
Tr[(— (h 2 /2)A — (a + V + ). 1^)070]. More precisely, since V < C (|x| _1 + 1), we have 



Tr[(-(/i 2 /2)A - (a + V+). 1^)070] > Tr[(-(/i 2 /4) A - C |x| _1 + 1)070] 

+ Tr[(-(/i 2 /4)A-(Q ) + a + l).ln 
> -Ch- 2 \\nh\-C(a + l) 2 h- 2 \n\. 



□ 



Proof of Lemma[]M 1- The lower bound follows from the radial symmetry of p and the fact 
that A(|x| _1 ) = |x| -3 > pointwise for all x 7^ 0. 
In fact, since p is radially symmetric we can write 
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where dz\ and dz 2 are normalized Lebesgue measure on the circles S x := {z G M 2 : \z\ = \x\} 
and S y := {z E R 2 : \z\ = \y\}. 

If \x\ > \y\ then using the subharmonic property of the mapping z i— >■ |x — in the 
open set {z e I 2 : |z| < we get 



^2 



Together with the similar inequality for \y\ > \x\, we obtain the desired lower bound on 



2. Because p(rr) < (27r|a;|) 1 and f p = \, for any k > 1, 
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Optimizing the latter estimate over k > 1 yields the first upper bound on p * | . 

3. We now prove the second upper bound on (p * | . | _1 )(a;) for \x\ large. We start by 
decomposing R 2 into three subsets 
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Fix e > small. For \x\ large enough, 
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Moreover, since p(y) < (27r|y|) 1 , 
(All) 
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Next, using the polar integral in 

Q 3 G{y eM 2 : 3\x\/2 > \y\ > \x\/2, \\x\ - \y\\ > \x\' 2 } 
we have, with notation s := min{|x|, r}/max{|x|, r}, 

P{V) dy< [ I dOdr. (38) 



\x-y\ J J max{|x|,r}v'l + s 2 - 2scos(6») 

s«3 3|x|/2>r>|x|/2,|r-|x||>|a;|- 2 

The singularity of the integral w.r.t. 9 (at s — > 1 _ ) is controlled by the following technical 
lemma (we shall prove later). 

Lemma 24 (Upper bound on elliptic integral). There exists a finite constant C > such 
that 

/dO 
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Using fl38l) and Lemma [241 we get, for \x\ large enough, 



/^L d ,< Cl Mai f p(y)iy (39 ) 



«3 3|:r|/2>|y|>|x|/2 



for some universal constant C7j. 

Putting (|36|) . (1371) and (1391) together, we conclude that for any e > there exists R = 
R(e,p) such that for any |x| > R, 

Fl Fl J 

3\x\/2>\y\>\x\/2 

□ 

For completeness we provide the proof of the upper bound on the elliptic integral. 
Proof of Lemma We just need to consider the singularity when s — > 1~. Write 

1 + s 2 - 2s cos(2fl) = (1 + s) 2 - 2s(l + cos(0)) = (1 + s) 2 - 4s cos 2 (#/2). 

Denoting k 2 = 4s/ (1 + s) 2 and making a change of variable (6 i— > 7r — 2$), we need to show 
that 

tt/2 1 

tf(Jfe) := / - M = [ dt = < C\ ln(l - k)\ 

J ^l-k 2 sm 2 (6) J ^(l-t 2 )(l-k 2 t 2 ) ~ 
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when k — > 1 . This upper bound follows from the identity 

dt 1 / 1 + 

7(1-^(1 -i*y ~ 7^ 11 1 1 - Vk 



□ 



Remark. The function if(fc) is the complete elliptic integral of the first kind. Its asymptotic 
behavior at k — > 1 _ is well known. It is (see [Tj, eq. (17.3.26), p. 591) 

K(k) = ^\ln(l-k)\ + hn(2) + o(l) k ^-. 

Proof of Lemma\2]\ Recall that we are working on the region 2A > \x\ > r ^> s. We start 
with the triangle inequality 

\V-V*gl\<\\^-\A-Ugl\ + \pY*\^-pY*\.\-^gl\. (40) 

(If V(x) = — 1 then the term invloved p^ F disappears.) 
When \x\ > r ^ s > \y\ using 

— | a; — 2/ 1 1 1 — Cs\x\~ 2 

one has 

li • r 1 - 1 • r 1 * g 2 s \ (x) < f i m- 1 -\ x - yr 1 ! g 2 Md y < c s \x\- 2 . (4i) 

Moreover, 

IpI f * i • r 1 - p7 * i • r 1 * 9 2 s \ w < || pI f (x - y ) i br 1 - 1* - ^r 1 ! ^)^- (42) 

We divide the integral into two domains. If \y\ > r/2 then using 

\\y\- 1 - \y - z\-'\ ^Cs^KCsr- 1 ^- 1 
and (pJ F * I . | _1 )(x) < C(|x| -1 / 2 + 1) (see Lemma IT2"|) we obtain 

pY(x - y) I \y\~ l -\y- z\~ x \ fi(z)dydz < Csr^ip™ * | . ^(x) < Osr^- 1 / 2 . (43) 
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If \y\ < r/2 then using 
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Replacing ( j45J) and (144 p into 042 p we arrive at 

|/°I F * I • l _1 — P\ F * I • l _1 * 9s\ ( x ) — C(sr~ 1 \x\~ 1 ^ 2 + r|x| _1 + r) when \x\ > r. 
From the latter inequality and (|4"T|) we can deduce from (|4"0"|) that 

\V * g 2 — V\(x) < C(s\x\~ 2 + sr~ 1 \x\~ 1 / 2 + r|x| _1 + r) when \x\ > r. 

Taking the square integral of the previous inequality over {r < \x\ < 2A} we get (with 
A = \ \nh\, r = h 1/2 , s = h 2/3 ) 



J T' 



\V-V* g 2 s \ 2 (x)dx < C(s 2 r- 2 A + r\ ln(A/r)| + r 2 A 2 ) < Ch 1/4 . 



r<\x\<2A 

□ 
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